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Abstract 


The  hypothesis  that  the  fouirth  moments  of  the  field  amplitude  distribution 
are  related  to  the  second  moments  as  in  a  normal  distribution  is  examined  to  determine 
whether  it  is  consistent  with  the  equations  of  motion  for  stationary,  hcwiogeneous, 
and  incompressible  hydroraagnetic  turbulence.  An  application  of  this  hypothesis 
to  the  two-point,  two-time  amplitude  distribution  of  the  velocity  field  and 
magnetic  field  leads  to  a  gross  violation  of  energy  conservation  at  high  Reynolds 
nimber.  This  result  holds  also  in  the  non-magnetic  case.  An  estimate  of  the 
magnitude  of  the  violation  yields  a  rate  of  energy  generation  of  order  E.  »  E,/C- j 
where  E,  is  the  energy  contained  in  the  inertial  range  and'J^  is  a  period 
characteristic  of  modes  in  the  inertial  range.  The  results  obtained  are  discussed 
in  relation  to  a  recent  theory  of  turbulence  formulated  by  Chandrasekhar. 
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1.  Introduction 

A  number  of  workers  in  turbulence  theory  have  employed  the  hypothesis  that 
certain  fourth  moments  of  the  statistical  distriVmtion  of  the  field  amplitudes  are 
related  to  second  moments  as  in  a  normal  distribution.  This  hypothesis  app^srentlv  . 
was  first  employed  by  Millionshtchikovl-  ^  .    Later,  it  was  applied  to 
the  two-point,  simultaneous  velocity  distribution  by  Heisenberg'-  ■',  Obukhoff  I-  -*, 
Batchelori'  -'j  Chandrasekhar'-  ■',  and  others,  in  order  to  estimate  pressure  fluctuations 
in  turbulent  flows,  Heisenberg'-  ^   used  an  implicit  application  of  the  hypothesis 
to  the  two-time  velocity  distribution  in  an  attempt  to  develop  equations  determining 
the  spectmjm  of  high  Reynolds  number  turbulence,  and,  more  recently,  a  theory  of 
turbulence  dynamics  employing  this  form  of  the  hypothesis  has  been  fonnulated  by 
Chandrasekhar '-  -'  and  extended  to  hydromagnetic  turbulence  ^-^ .  Proudman  and  Reid^  -^ 
had  previously  applied  the  hypothesis  of  normal  relationship  to  the  three-point, 
simultaneous  velocity  distribution  in  a  treatment  of  the  decay  of  homogeneous 
turbulence.  Some  experimental  evidence  with  regard  to  these  hypotheses  is  discussed 
by  Uberoi  •-  J  and  Batchelor  ^  -^ . 

In  the  present  paper,  the  consistency  of  the  hypothesis  of  nonnal  relationship 
of  fourth  and  second  moments  is  examined  for  the  esse  of  stationary,  homogeneous, 
incompressible  hydromagnetic  turbulence.  It  is  concluded  that  this  hypothesis 
applied  to  the  two-time,  two-point  distribution  of  the  field  amplitudes  leads  to 
gross  violation  of  energy  conservation  at  high  Reynolds  numbers.  The  conclusions 
reached  are  discussed  in  relation  to  the  dynamical  theory  of  Chandrasekhar,  mentioned 
above. 

2o  The  equations  of  motion 

Consider  an  incompressible  fluid  of  uniform  density  p,  kinematic  viscosity  v, 
permeability  n,  and  conductivity  a  ,  which  is  subjected  to  an  external  force  field 
F(x,t),  Under  the  related  assumptions  of  nonrelativistic  motion,  negligible  dis- 
placement current,  and  high  conductivity'-  J,  one  may  derive  from  Maxwell's  equations. 
Ohm's  law,  and  conservation  of  the  total  stress-energy  tensor  the  divergence  conditions 

(2.1)  V*B  =  7-u  ■  0 

and  the  equations  of  motion 

(2«2)  B  -  V  X  (iucB)  -  (Untitr)""'-  c^  V^  b*  -  0 

(2»»  i   +  (u»7)u  -  (Unup)"-'-  (Vxb)  x  B  +  p"-""  Vp  -  vV^  u  -  F, 

where  ulx,t)  is  the  velocity,  B(x,t)  is  the  magnetic  induction,  and  p(x,t)  is  the 
pressure. 
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Equations  (2.2),  (2.3)  may  be  solved  for  p  and  rewritten  in  the  tensor  form 

(2.U)  w.  -  vV  w.  »  7.(u.w.  -  u.w.  ) 

^  1      1    J  1  J    0  1 

(2.5)  \   -  vV^  u.  =  Pij(^)[^|(V/  "  ""jV  *  ^J  ' 
where 

(2.6)  w  -  (Unup)"^/^B,  V  -  (Un^o-)"-^  C^, 


and 


(2.7) 


^ j  ^'^  ■  ^io  -  '■'  'i':5         ['i  ■  ^/^  ^i 


The  projection  operator  P.  .(V)  applied  to  the  nonsolenoidal  part  of  a  vector  field 
gives  zero^  it  represents  the  action  of  the  pressure  forces.  Both  divergence  condi- 
tions (2.1)  are  preserved  by  the  equations  of  motion. 
The  total  energy  of  the  system  is 

(2.8)  Energy  =«  5-  pj  (  u^  +  w^Jd-'x. 


^  p/{  u^  .  w2]< 


If  the  flow  is  of  finite  extent,  or  obeys  cyclic  boundary  conditions,  one  may 
verify  by  partial  integration  the  relation 

(2.9)  /K'j^Vj  -  Vi^  "  ^^j^'^!£^^j"jt  -  "j^jL^]^^^  ■  °' 

which  demonstrates  the  conservative  nature  of  the  nonlinear  interaction. 

It  is  convenient  in  what  follows  to  express  the  field  amplitudes  in  Fourier 
series,  employing  the  artifice  of  cyclic  boundary  conditions  on  the  faces  of  a  cube. 
An  equivalent  analysis  could  be  carried  out  in  terms  of  Fourier  integrals,  but  at 
the  expense  of  more  cumbersome  notation.  The  equations  of  motion  for  the  Fourier 
coefficients  associated  with  wave  vector  k  are 

(2,30)  k^w^(ft  -  k^u.(k)  -  0 
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(2.11)  v^(k)   +  vk^  w^(k)  -  ikj^jf,    LciT-  k«)w  (iT')  -  u.(k  -  l^')w^(k')] 

(2.12)  u^(k)  +  vk^u^(k)  -  ik.Pj^j(k)  Z^-^  ,  rw^(k  -  k')w£(t')  -  Uj  (k  -  k')u£(k') 


*P,j(k)F.(k), 


and  the  energy  per  unit  mass  is 


1       r ""  -^     -i 

(2.13)  ^  ^  2  H  rL^i^^K^^^   '  ^(k)w.  (k) 


Let  the  linearly  independent  real  and  imaginary  parts  of  the  independent  vector 
components  of  all  the  Fourier  coefficients  u(k),  w(k)  be  arranged  in  some  one- 
dimensionai  serial  order  and  be  denoted  by  q  //5",  a  ■  1,2....  For  each  allowed 
wave  vector  k  there  are  two  independent  vector  components  of  each  of  the  Fourier 
coefficients  of  the  real,  solenoidal  vector  fields  u  and  w;  thus,  there  are  eight 
q's  for  each  allowed  pair  k,  -  k. 

The  equations  of  xotion  of  the  q's  may  be  written 

(2olU)  q^  ^  v^q,  -  E  P,y^Py^pS-  *  ^a    [^apr  "  W]' 

where  the  coefficients  A  „  and  the  damping  factors  v  contain  the  k's  implicitly 
and  the  f  are  the  suitably  ordered  and  normalized  real  and  imaginary  parts  of  the 
independent  vector  components  of  the  solenoidal  force  field  P.  .(7)F..  All  the 
quantities  appearing  in  (2.1U)  are  real,  and  the  v  are  all  positive.   (No  sum- 
mation convention  is  employed  for  the  Greek  indices.) 
The  total  energy  per  unit  mass  is  given  by 


(2.15)  E-|Ea<\)'' 


and  equation  (2.9),  which  expresses  conservation  of  energy  by  the  nonlinear  inter- 
action, implies 


(2-i«       H  <,,p,AprVBS- 


0. 


-  u- 

This  relation  must  hoid  for  all  values  of  the  independent  variables  q.  Taking  all 
the  q's  to  be  zero  except  q  ,  q^,  q  ,  one  finds 

(2.17)        A^^^  ^^^^^  V^.x'°• 

If  the  three-mode  interaction  of  q  ,  q..,  q^  is  defined  by  those  nonlinear 
terms  in  the  equations  of  motion  of  these  three  valuables  which  contain  no  other 
q's,  it  is  apparent  from  (2.1U),  (2.17)  that  this  elementary  interaction  contributes 
a  vanishing  part  to  the  time-rate  of  change  of  E  and  is  therefore  individually 
conservative. 

The  selection  rules  to  determine  ^ich  A's  are  nonvanishing  are  fixed  by  the 
orthogonality  conditions  (2.1)  and  by  which  wave  vector  pairs  contribute  to  the 
convolution  sums  in  equations  (2.11),  (2.12).  It  will  be  assumed  that  there  is 
no  mean  flow  and  no  uniform  component  of  magnetic  field.  Then  the  Fourier  coeffi- 
cients u(0),  w(0)  vanish,  and  the  orthogonality  conditions  permit  only  pairs  of 
Fourier  coefficients  such  that  k  and  k  -  k'  are  nonparallel  to  contribute  to  the 
convolution  sums.  In  this  case,  the  selection  rules  on  the  A's  include  the  condition 

(2ol8)  A  ,   *  0,    unless  m-,X,  cr  are  all  different. 

If  the  Fourier  coefficient  F(0)  vanishes,  the  conditions  u(0),  w(0)  ■  0  are  preserved 
by  the  equations  of  motion. 

The  various  equations  of  motion  (2.U),  (2.5),  (2.11),  (2.12),  (2.1U)  may 
be  integrated  with  respect  to  time  to  yield  integral  equations.  From  (2.5)  and 
(2.1U)  one  obtains,  respectively, 

(2.19)  u^  «  [   P^j(7)exp[v(t-f)V^][7^(wjWj'-u^u^)  +FJdt' 

(2.20)  q^  .  j   exp[-v(t-t . )]  [H  A^^^q-q^  *   f 'Idt', 

Where  unprimed  variables  have  argument  t,  primed  variables  have  argument  t',  and 
the  boundary  condition  of  finiteness  at  t  »  -  oo  has  been  applied. 

3.  The  statistical  equations  of  homogeneous  stationary  turbulence 

The  statistical  theory  of  turbulence  is  concerned  with  distributions  of 
solutions  of  the  equations  of  motion  in  a  function  space  of  space-time  functions. 
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It  is  usually  presumed,  on  physical  grounds,  that  the  finite-order  moments  of  the 
distribution  and  their  finite-order  derivatives  exist.  Stationary,  homogeneous, 
or  isotropic  distributions  are  defined  by  the  invariance  of  the  moments  under  time 
displacement,  space  displacement,  or  proper  rotations,  respectively. 

In  the  stationary,  homogeneous  case  to  be  considered  here,  one  may  define  the 
moments 


f^    C.j(y,r)  =^u^(x,t)u.(x',t')> 
fi^j(y,r)  =^u^(x,t)u^(xSf)^ 


(3.1) 


V^ 


^ijk^y^'"^^  "  <u^(x,t)u^(x,t)uj^(x',t')^ 
\^yiifl:^y'^)   =^i(it)uj(x,t)u^(^t')u^(x,t')/' 
Gij(7,r)  =<F.(x,t)Uj(x;t.)> 
%  jj^ (y,'Z^)  -  <(u^  (x,t  )uj  (x,  t  )F^  (xjt '  )> 
[y  -  X-  X',  er-  t  -  f]  , 


which  belong  to  the  two-point,   two-time  distribution  of  u  and  F;   similar  moments 
involving  w  may  also  be  defined.     It  follows  from  the  definitions  that  Z.J^yZ)  » 
C..(-y,-'t),   with  analogous  relations  for  other  moments. 

Multiplying  equation     (2.5)  for  u.  (x,t)  by  u,  (x,'t   )  and  averaging  over  the 
distribution,   one  obtains,   for  the  special  case  w  ■  0,  the  statistical  equation  of 
motion 


(3.2) 


^ik-^'\k"Pij('>r- Vjxk^^jki' 


,-»»i  I  • 


where  V  operates  with  respect  to  y.  If  u,  (x, t  )  is  expressed  according  to  the 
integral  equation  (2.19)  before  the  averaging,  one  can  obtain,  instead,  the  integro- 
differential  equation 
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Cilj(y,0)  -  vW^C^j^(y,0)  =  (  P   (V)V^exp(vtV^) 

(3.3) 

^  L-  Pk:n(^)'n^j/mn(y»-^>  ^  H^,(y,-r)]dr.  P.  .  (V)G^^(5^,0), 

in  which  the  fourth  moments  £.«   appear  instead  of  the  third  moments  D...  . 

jjfran  "  j^K 

The  moments 


V^^'<W"^^ 


(3.U) 


describe  the  two-time  q-distribution.  The  definitions  require  R  .  ("2')  ■  R^  {-T) 
with  analogous  relations  for  the  other  moments.  In  analogy  to  equations 
(3.2),  (3.3),  one  may  derive  the  statistical  equations  of  motion 


CD 

(3.6)  R  (0)  +  V  R  (0)=  5^   A^A-   I    e  ^'^'^T.  ^    {'l;)dt, 

valid  for  modes  such  that  f  "0. 

The  quantity  R  (0)  ■  ^q  q  \  represents  the  average  rate  of  increase  of  the 
energy  of  the  mode  q  and  must  vanish  in  a  stationary  distribution.  It  is  apparent 
from  (3.5)  that  the  contribution  to^q  q"^from  the  elementary  three-mode  inter- 
action of  distinct  modes  q  ,  q^,  q,^  is 

(3.7)  J\   ^p^  =  2A^pY  Sp^(O)  .  2A^p^<q^qpq^>  . 

The  indi\'idual  conservation  property  of  the  three-mode  interactions  expressed  by 
equations  (2,1?)  has  the  consequence 
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(3.B)  ^apY*  -^PYa  "    ^YcS  =  °- 

For  fields  u  and  w  obeying  cyclic  boundary  conditions  on  the  faces  of  a 
cube  of  side  L,  the  requirement  of  invariance  of  the  statistical  distribution 
under  spatial  displacement  modulo  L  implies  that  the  Fourier  coefficients  associated 
with  different  wave  vector  pairs  k,  -k  must  have  zero  correlation  and  that  the  real 
parts  of  the  coefficients  for  a  given  wave -number  pair  must  be  uncorrelated  with 
the  imaginary  parts.  In  addition,  it  is  always  possible,  in  this  case,  to  choose 
principal  axes  for  each  coefficient  so  that  the  two  independent  vector  components 
of  u(k)  are  uncorrelated  and  the  two  independent  components  of  w(k)  are  uncorrelated  . 
The  equations  of  motion  (2.U),  (2.5^)  are  invariant  under  sign  reversal  of  w.  This 
ensures  the  existence  of  solution  distributions  in  which  u  and  w  are  uncorrelated, 
as  has  been  noted  by  Chandrasekhar  •-  -^ .   Since  the  q's  are  the  independent  vector 
components  of  the  real  and  imaginairv  parts  of  the  Fourier  coefficients  u(k),  w(k), 
it  follows  that  for  such  distributions  the  requirement  of  hcanogeneity  and  the  choice 
of  principal  axes  described  yields  the  conditions 

(3.9)  R^p(tr)  -  6^pR^Jt). 

A  two-time  distribution  will  be  called  quasi -normal  if  for  every 
choice  of  fovir  quantities  a,  b,  c',  d'  (where  no  prime  denotes  time  t  and  a  prime 
denotes  time  t')  the  fourth  and  second  moments  are  related  by  the  equation 

(3.10)  <abc'd'>  =  <ab>  <c'd«>  +  ^c'>  <bd«)  +  ^d;>  <(bc '^  , 

which  holds  for  a  normal  distribution.     The  hypothesis  that  the  two-time  stationary 
distribution  of  the  q's  is  quasi-normal  implies 

If       A  -  y"a  q  ,    B  =  TB  q  ,   C     =  5~Y  ^'j    D     =  5~5  q'     are  any  four  linear 

combinations  of  the  q's  with  coefficients  constant  over  the  distribution, equation 
(3.n)  yields 

^— ^KC— ^i^"^^»^— —  ■  ■  ■■  - t^     •  ■!■  — .1  ■■■■11  II      ■     ■   —  —  —  ■  I— . ___^ . . . ■■  ^ 

For  isotropic  turbulence,  this  is  true  independently  of  the  choice  of  axes.  Symmetry 
conditions  on  the  correlation  and  spectrum  tensors  in  homogeneous  turbulence  are 
discussed  in  [lo].  Chapter  III. 
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(3.12)  {abc'd^  -  ZVxV^2'<VVa^r> 

-    ^b^<^'d'^  +  /1c'\/bd'^  ♦  <ad'^<(bc'^     . 

Thus,   since  the  q«is  are  the  real  and  imaginain'  parts  of  all  the  Fourier  coefficients 
of  u  and  w,   equation  (3.11)  implies  that  the  four-point,  two-time  joint  distribution 
of  u  and  w  is  quasi-normal. 

The  hypothesis  that  the  two-point,   two-time  joint  distribution  of  u  and  w  is 
quasi-normal  is  evidently  considerably  less  restrictive  than  equation   (3.11). 
However,   it  is  important,   in  what  follows,  to  note    that  if  any  two-point,   two- 
time  fourth  moment  of  u  and  w  is  reduced  to   second  moments  by  the  two-point  quasi- 

normality  hypothesis,    and  then  expressed  in  terms  of  the  R  „  by  Fourier  analysis, 

ap 

the  result  must  be  identical  with  that  which  would  be  obtained  by  applying  (3.11) 
without  restriction.  This  may  be  seen  by  considering  u  and  w  as  sums  over  the  q's 
and  employing  (3.12)  in  reverse.  The  hypothesis  that  the  two-time,  two-point  joint 
distribution  of  u  and  w  is  quasi-normal  implies 

(3.13)        ^ijki^y*'^^  "  C^^(0,0)Cj^^(0,0)  ♦  C^j^(y,^)Cj^(y,^) 
with  similar  conditions  for  moments  involving  w. 


U.  Consequences  of  the  quasi-normality  hypothesis 

It  may  be  seen  very  simply  that  the  quasi-normality  hypothesis  (3.11)  is  incon- 
sistent with  the  conservation  properties  of  the  nonlinear  interaction.  Consider  any 
three  modes  q  ,  q^,  q     such  that  f  ,  f„,  f  vanish.  They  could  be  magnetic  modes 
or  velocity  modes  in  the  inertial  or  dissipation  ranges.  Substituting  the  integral 
expression  (2.20)  for  q  in  (3.7),  one  finds  for  the  mean  rate  of  transfer  of 

energy  to  mode  q  through  the  three-mode  interaction  with 
a 


y  V 
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(U.l)  -/I  ^  -  2  A  -  y;   a  ,   f  e  "^  T  .„  ('D)d^ 

'o 


Y^  11,  X  anX  J 

'0 

If  the  hypothesis  (3.11)  is  asserted,  (U.l)  may  be  rewritten*^" 


where  the  relations  (2.18),  (3.9)  have  been  used. 

It  may  now  be  noted  that  the  integrals  on  the  right  side  of  (U.2)  are  all 
positive  quantities.  This  follows  from  the  fact  that  the  spectrum  of 

e       is  everywhere  positive  and  the  spectra  of  the  correlation  functions 

Roo('2')»  R  (^  cannot  be  negative  anywhere.  Since  the  spectrum  of  a  product  of 
pp      YY 

functions  with  non-negative  spectra  has  itself  a  non-negative  spectrum,  and  since 

the  integral  of  the  even  function  e      R  q('C)R  CC*  )  gives  the  value  at  the 

origin  of  its  Fovirier  transform,  none  of  the  integrals  can  be  negative.  Given 

V  ^  0  it  follows  further  that  none  of  the  integrals  can  be  zero,  if  R„q(C'), 
a  pp 

R (-C)  are  well  behaved.  This  may  be  seen  by  noting  that  R--(O),  R.^(0)  >  0, 

YY        .^j^i  PP      YY 

so  that  if  e       is  replaced  by  a  new  function  in  which  the  cusp  is  rounded 
off  arbitrarily  close  to  the  origin,  but  whose  spectrum  is  still  everywhere  non- 
negative,  the  values  of  the  integrals  are  reduced  but  must  still  be  non-negative. 
If  one  takes  v  "  0,  a  given  integral  can  vanish  only  if  the  corresponding  pair 
of  modes  have  no  frequencies  in  common  over  the  distribution,  that  is,  only  if 
there  are  almost  no  pairs  of  solutions  in  the  distribution  such  that  q^  of  one 
has  frequencies  in  common  with  q_^  of  the  other. 

Applying  the  same  procedure  and  argument  to  iL  „   and  _/T.  „,  we  obtain 

(U.3)  napY>°'  -^PYa>°*  ^YaP  =^  °* 

which  violates  the  detailed  conservation  condition  (3»8). 

Certain  consequences  of  the  two-point,  two-time  quasi-normality  hypothesis 
will  be  discussed  now.  Consideration  will  be  confined  to  the  nonmagnetic  case 
w  «  0.  This  results  in  appreciable  reduction  in  manipulations,  and  the  generaliza- 
tion to  the  magnetic  case  (with  appropriate  augmentation  of  equation  (3. 13)) is 
straightforward. 


■^The  use  of  integrated  equations  of  motion  to  relate  fourth  and  second  moments 
apparently  was  first  described  by  W.  HeisenbergL'^J  . 
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The  decomposition  of   (3.3)  into  Fourier  series  3rields  an  expression  for  the 
Fourier  coefficients  of  C,  (y,0)  in  terms  of  the  Fourier  coefficients  of  the 
tensors  3 ,  A  y  ^  t  7^ '     Consider  a  wave  vector  k  such  that  the  Fourier  coeffi- 
cient F(l^)  of  the  driving  force  vanishes.     The  Fourier  coefficients  of  ^  and  rr 
are  cross-power-spectra  of  f^(x,€^  with  functions  of  u(x,t)  and  must  vanish  for 
such  P^.     The  Fourier  coefficients  of  C^ij(y»0)  are   <Cu*(k)uj^(k)^   and  determine 
R     (0),  if  q     is  a  velocity  mode  belonging  to  wave  vector  k.     It  follows  that 
equation  (3.3)  permits  us  to  express  R     (0)  in  terms  of  Fourier  coefficients  of^ 
and  ^  alone  if  q     is  a  velocity  mode  associated  with  a  wave  number  for  which 
F(k)  -  0. 

If,   now,   the  two-point,   two-time  quasi-normality  hypothesis   (3.13)  is 
asserted,  the  tensor  ^  may  be  expressed  in  terms  of  the  tensor  o  ,   and,   in 
consequence,  the  Fourier  coefficients  of  ^  may  be  expressed  in  terms  of  sums 
of  products  of  Fourier  coefficients  of  ^   and,   in  turn,   expressed  in  terms  of  the 
quantities  R  „.     According  to  the  discussion  in  the  preceding  section,   this 
expression  in  terms  of  the  R  „  must  be  identical  with  that  obtained  by  applying 
the  stronger  quasi-normality  hypothesis   (3.11)  to  the  equations  of  motion. 
Applying   (3.11)  to   (3,6)  aixl  specializing  to  the  case  a  ■  Hj   one  finds  with  the 
aid  of  the  relations   (2.18),    (3.9), 


.00 


(U.U)  R^JO)  .  2  E  (A,p^)'    {       e"'-^Rpp(tr)R^(  t  )d  t  -v^R,,(0). 

The  quantity  R  (0)  ■  ^q  q  ^  is  the  average  rate  of  increase  of  the  energy  of 

the  mode  q  and  must  vanish  in  stationary  turbulence.  The  vanishing  of  the  right 
a 

side  of  (U.U)  is  the  condition  of  energy  balance  between  the  net  energy  input  due 
to  nonlinear  interaction  and  the  viscous  dissipation  for  a  mode  not  subjected  to 
external  driving  forces. 

It  is  not  immediately  apparent  that  (U.U)  is  impossible  to  satisfy,  but  the 
relation  may  easily  be  seen  to  indicate  behavior  at  variance  with  the  usual  physical 
picture  of  the  dynamics  of  the  energy  cascade  process  in  high  Reynolds  nximber 
turbulence.  Consider  a  mode  in  the  inertial  range.  It  is  accepted  that  the  dynamics 
of  such  modes  are  almost  entirely  controlled  by  the  Rejmolds  stresses  and  pressure 
forces  and  that  the  direct  effect  of  viscosity  in  this  range  is  negligible.  If 
one  now  imagines  reducing  the  damping  factor  v  while  leaving  undisturbed  the 
damping  factors  of  the  other  modes,  the  effect  on  the  excitation  of  q  and  on  the 
dynamics  of  the  other  modes  should  be  negligible.  According  to  (U.U),  however,  if 


-  11  - 

• 

R  (O)  is  to  remain  zero  the  effect  of  reducing  v  by,  say,  a  factor  of  10  must 
be  profound.  Either  the  mean  excitation  R  (0)  of  the  mode  q  wauld  have  to 
increase  by  a  large  factor,  or  the  positive  integrals  on  the  right  would  have 
to  change  significantly. 

In  any  individual  solution  representing  a  turbulent  flow  at  high  Reynolds 
number,  the  time  variation  of  a  mode  q  in  the  inertial  range  should  be  exceedingly 
COTiplicated  and  its  frequency  spectrum  should  be  broad.  In  a  flow  not  subjected 
to  random  external  perturbations  it  is  perhaps  conceivable  that  the  frequency 
spectra  of  modes  in  the  inertial  range  could  have  some  rapidly  oscillating  structure, 
but  in  an  actual  flow  such  oscillations  surely  would  be  smeared  out  as  a  result 
of  perturbations.  It  follows  that  in  a  distribution  representing  the  presumed 
ergodic  behavior  of  high  Reynolds  number  turbulence,  the  averages  Rp^C'C),  R     {t) 
should  be  smooth  well-behaved  functions  of  ^  .  Further,  it  seans  likely  that  thf'Se 
functions  should  not  exhibit  significant  regions  of  negativity.  This  connnent  is 
based  on  the  observation  that  the  forces  A.  ,q  q^  exerted  on  mode  q-  through  the 

\i[iK   [1.   K  p 

three-mode  interactions  are  quadratic  in  the  q's,  so  that  their  frequency  spectra, 

and  consequently  the  frequency  spectrum  of  q„  should  contain  strong  zero  frequency 

^oo         ^ 
components;  hence  the  moments  |   R  ("C)  d  C  should  not  differ  greatly  in  order 

of  magnitude  from  CqRqq(O),  where  2'.  is  some  characteristic  period  for  the  mode. 

p  PP  P 

If  R  (^),  Roo('C'),  R  ('D  in  (U.h)  are  such  smooth  functions,  it  is  difficult 
aa      PP      Tf 

to  see  how  they  could  adjust  to  a  large  reduction  of  v  without  significant  change 
of  the  excitation  levels  or  the  characteristic  mode  periods.  In  the  next  section 
an  estimate  will  be  made  of  the  magnitude  of  the  sum  over  integrals  in  equation 
(U»U)  for  high  Reynolds  nioraber. 

5.  Discussion  of  a  nondissipative  system 

An  insight  into  the  relation  of  the  quasi-normality  hypotheses  to  conservation 
of  energy  by  the  nonlinear  interaction  is  provided  by  the  study  of  an  idealized 
flow  system  in  which  both  dissipation  and  external  forces  are  assumed  to  vanish.  As 
in  the  last  section,  explicit  consideration  will  be  confined  to  the  nonmagnetic  case, 
w  "  0,  the  extension  to  the  more  general  case  being  easy  but  more  cumbersome. 

If  V,  F  and  w  are  taken  equal  to  zero,  it  is  well  known  that  the  solutions  of 
(2.5)  develop  surfaces  of  discontinuity  which  correspond  physically  to  shock  fronts. 
It  is  not  possible,  for  this  reason,  to  speak  meaningfully  of  nontrivial  stationary 
distributions  of  solutions  of  the  inviscid  equations.  To  avoid  such  behavior,  the 
altered  equation  of  motion 
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(5.1)  u.  -  -»^(K"272)p^j(V)V^  (u.u^) 

will  be  adopted.  The  operator  n  ,  defined  by 

(5.2)  )^(S)  -  1,  |S!  <  1,  K|(S)  »  0,  |SI  >  1, 

prevents  the  development  of  infinite  gradients.  It  may  be  verified  by  application 
of  Green's  theorem  that  the  property  of  conservation  of  the  energy  expression 

(2.3)  is  not  altered  with  arbitrary  choice  of  the  operator  y\. 

In  terms  of  the  Fourier  representation,  h  provides  a  cut-off  in  wave  number, 
suppressing  completely  all  modes  of  wave  vector  |kl  >  K.  It  may  be  regarded 
"physically"  as  representing  an  infinite  viscosity  for  these  modes  .  The  equations 
of  motion  for  the  finite  number  of  modes  which  remain  are 

(5.3)  \-ZVy^pS-   t£-l»  2,  ...¥), 

where  N  is  the  total  number  of  q's.  The  A's  may  be  taken  as  identical  with  those 
in  (2.1b)  if  all  A's  referring  to  wave  numbers  |k  |  >  K  and  to  magnetic  degrees  of 
freedom  are  removed.  The  detailed  conservation  property  (2.17)  provides  another 
demonstration  that  the  operator  ^  does  not  destroy  energy  conservation. 

It  is  evident  from  equations  (5»3),  (2.18)  that  3q  /dq  =  0,  so  that  the  new 
system  obeys  the  Liouville  equation 

N   aq 

(5.'')         X  ST  "  °- 

a«l   ^a 

An  immediate  consequence  is  that  any  distribution  of  simultaneous  amplitudes  which 
depends  only  on  the  constant  of  notion  E"^^  (q)  is  time  invariant  under  the 
equations  of  motion;  for  such  distributions  there  is  equipartition  of  energy  among 
the  modes.  In  particular,  the  Maxwell  distribution  of  simultaneous  amplitudes 

(5.5)  Z  exp  [-P  Xi^^a^^]' 

where  Z  is  a  normalization  constant  and  p  is  related  to  the  mean  energy  per  mode. 


ft 

Application  of  a  similar  cut-off  to  an  inviscid,   adiabatically  compressible  fluid 
is  discussed  in  [l2] . 
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is  time  invariant. 

If  the  hypothesis  that  the  two-point,  two-time  distribution  of  u  is  quasi- 
normal  is  now  asserted,  the  procedures  of  the  last  section  lead  to  the  equation 

OD 

(5.6)  E  .^  q^q^>  -  27:  %/    (        ^p^C  ^  )R^(  t  )d  T  . 

Energy  conservation  requires  that  £  vanish.  As  discussed  before,  this  would  be 

possible  only  if,  given  q  and  q.^  such  that  A  .  /  0,  (almost)  no  functions  q 

P      I  QpY  p 

associated  with  any  given  solutions  in  the  distribution  should  have  any  frequencies 

in  common  with  the  functions  q^^  associated  with  (almost)  any  other  solutions  in 

the  distribution.  This  would  be  absurd  . 

It  is  possible  to  obtain  a  crude  estimate  of  E  in  (5.6)  if  one  adopts  the 

simultaneous  amplitude  distribution  (5«5)  and  assvimes  that  the  distribution  averages 

R  (TT)  are  smooth,  essentially  positive  functions  for  N  large.  Assume  the  modes 
ao. 

are  coupled  in  such  a  way  that  all  the  characteristic  mode  periods  f    are  of  the 
same  order  of  magnitude  ■? •  Let  the  mean  energy  of  the  equally  excited  modes  be 
K-R  (0)  •  s-R  .  Let  some  average  value  of  the  non-vanishing  coefficients  A  ^  be 


A  and  let  the  average  number  of  non-vanishing  A  „„  for  fixed  a  be  M.  Then  an  estimate 


(5.7)         f-'~  R-^  <4,5,>  -  R-J^Z:  ^  *apA.X^prux'°' 

where  the  relations  (2.18)  and  the  rigorous  independence  of  different  simultaneous 
q's  have  been  used.  Inserting  this  estimate  of  t    into  equation  (5.6),  and  estimating 


r 


R„„(r)R  (r)  dr  ~-  i-R^r  ,  one  finds 
PP     YY  2  ' 


(5.8)  t  E  ~  t^  NMA^  R^  -  I  NR  -  E 

for  the  change  of  energy  in  one  characteristic  period.  Thus  the  violation  of  energy 

Consider,  for  example,  the  system  of  just  three  modes  with  A^^.  ■  1,  A.,  •  -i, 
A,-_  '•  -s-  .   Its  equations  of  motion  are  invariant  under  exchange  of  q_  and  q.. 
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conversation  is  gross  . 

The  invlscid  system  just  described  may  be  considered  a  model  of  the  equilib- 
rium statistical  mechanics  of  the  inertial  range  of  high  Reynolds  number  turbulence. 
It  is  hard  to  see  how  the  transition  from  equilibrium  to  non-equilibrium  statistical 
mechanics  (non-equipartition),  or  a  more  careful  accounting  of  the  detailed  structure 
of  the  A's,  could  alter  the  estimate  (5.8)  by  a  very  large  factor.  The  result 
obtained  then  suggests  rather  strongly  that  if  the  hypothesis  of  quasi-normal 
two-point,  two-time  distribution  is  applied  to  actual  stationary  turbulence,  the 
net  transfer  of  energy  to  modes  in  the  inertial  range  indicated  by  the  sum  over 
integrals  in  equation  (U.U)  represents  a  total  rate  of  input  E.  -^  E./  f ^,  where 
E.  is  the  energy  contained  in  the  inertial  rango  and  f.  is  the  order  of  some 
characteristic  time  in  this  range.  The  viscous  dissipation  in  this  range  is  known 
to  be  very  much  less  for  high  Reynolds  nimiber. 

It  may  be  of  interest,  in  concluding  this  section,  to  outline  directly  in 
terms  of  the  field  u(x,t)  an  equivalent  of  the  argument  leading  to  (5.6).  The 
analog  of  equation  (3.3)  for  the  altered  non-dissipative  system  is 

00 

(5.9)  C,,(y,0)  -  I   [  ->!  2(^-2  V2)P.  .(^P^(V)y^  vJe.^  mn^^^'"^^^^- 

o 

Contracting  with  respect  to  i ,k,  and  setting  y  =  0,  one  obtains  on  the  left  side 
<\U.u.^  ,  which  is  proportional  to  the  mean  rate  of  increase  of  energy.  Under 
the  two  point,  two-time  quasi-normality  hypothesis  (3.13)  one  then  obtains 


(5.10)         < 


^i"i>'  I 

c 

[ 


00 


-  [=™j«.^'<=„i(--5^.'^  ♦  <^ni^-y'^K,^-y-^l .  ^^^' 


0 

where  the  identity  P. .P.  «  P.  has  been  used.  The  proof  that  the  right  side  of 
•^  ij  im  -  Jm  ^  * 

(5.10)  is  non-negative  may  be  shown  to  follow  from  the  fact  that  products  of  tensor 

correlation  fxinctions  are  themselves  correlation  f\mctions  and  from  the  fact  that 

2-2  2 

*J  (K   V  )P.  (y),  -7/V  are  positive  definite  Hermitian  tensor  operators.  Let  us 


00 


It  may  be  noted  that  even  if  the  moments  (  ^qqC^   ^"^  were  to  vanish,  the  quantities 

"  2  — 

R„ Q  ('r)R^^(t')dT  should  not  be  less  than  R  'T  by  a  large  factor,  since  the  mean  value 
pp    Tf 

O 

of  |RQQ('r)|  in  the  negative  region  should  certainly  be  less  than  the  mean  value  in 

pp 
the  positive  region  nearer  the  origin. 


]■ 
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denote  the  whole  integrand  in  (5»10)  by  Y(T)j  it  is  apparent  that  by  squaring 
{^»1 )  and  averaging  over  the  distribution  we  obtain 

(5.11)  {u^u.)  -  Y(0). 

This,  of  course,  cannot  vanish,  except  trivally,  so  that  if  the  dependence  of  YiT") 

on  '^  is  of  the  sort  supposed  above  for  RqoC'S'),  the  right  side  of  (5.10)  does  not 

pp 
vanish.  The  estimate  (5.8)  is  readily  corroborated  from  (5elO)  and  (5.11). 


6.  Conclusions 

The  investigation  carried  out  in  Sections  h   and  5  suggests  strongly  that  the 
two-time,  two-point  quasi -normality  hypothesis,  applied  to  high  Reynolds  number 
turbulence,  leads  to  a  gross  violation  of  energy  conservation.  Evidently,  normal 
relation  of  second  and  foiorth  moments  of  the  two-time,  two-point  distribution  does 
not  permit  those  distribution-averaged  phase  relations  among  coupled  triads  of  modes 
which  are  required  to  maintain  mean  energy  balance  in  all  the  elementary  three-mode 
interactions  . 

The  consequences  would  appear  to  be  serious  with  respect  to  the  foundations 
of  Chandrasekhar's  recent  turbulence  theory'-  -l*"-']^  which  is  based  on  this  hypothesis. 
However,  this  does  not  in  itself  imply  that  the  results  of  Chandrasekhar's  theory 
are  necessarily  badly  in  error.  The  consequences  of  the  quasi-normality  hypothesis 
depend  very  strongly  on  the  particular  way  in  tdiich  it  is  employed.  If  in  (3.7) 
one  substitutes  an  integral  expression  of  the  form  (2.20)  for  q.  or  q,^  instead  of 
q  it  is  readily  verified  that  an  essentially  negative  result  for  ./V  .  is  obtained. 
Thus,  a  suitably  symmetric  application  of  the  quasi-normality  hypothesis  couid 
result  in  expressions  for  the  energy  transfer  which  are  less  grossly  in  error  than 
equation  ih»k)     .  A  priori,  it  is  not  implausible  that  Chemdrasekhar 's  treatment 


The  set  of  moments  defining  any  realizable  statistical  distribution  must  obey  a  set 
of  generalized  Schwarz  inequalities.  In  the  case  at  hand,  if  second  and  third  moments 
of  the  two-time  distribution  of  the  q's  are  related  by  the  equations  of  motion,  it  is 
impossible,  under  the  quasi-normality  hypothesis,  to  satisfy  non-trivially  certain  of 
these  inequalities  relating  second,  third,  and  fourth  moments.  In  effect,  the  corre- 
lations between  q^  and  q^q^  and  between  \\(^aer» ^a^^'''  °'  ^'^  ^  ^''^^  required  by 
the  equations  of  motion  induce  a  correlation  between  q^q.^  and  q  q  so  that  the  moment 
^q^q^^q'q^^  does  not  vanish,  as  it  would  in  a  quasi-normal  distribution. 

a* 
Writing,  for  the  idealized  dissipationless  system, 

and  substituting  integral  expressions  for  q  ,  q.,  q^  respectively  in  the  three 
terms,  equation  (3oll)  gives  fl^p^  *  0  if  R^^(r)  -  Rpp(^)  -  ^w^^)- 
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in  terms  of  the  defining  scalars  of  isotropic  turbulence  accomplishes  this.  It  may 
not  be  concluded,  therefore,  from  the  evidence  presented,  that  Chandrasekhar • s 
theory  necessarily  leads  to  seriously  incorrect  physical  results.  What  does  seem 
evident  is  that  the  two-point,  two-time  quasi-normality  hypothesis,  which,  to  start, 
is  highly  arbitrary,  must  be  considered  a  most  suspect  theoretical  foundation  if 
only  a  particular  choice  of  the  method  of  its  application  givos  the  possibility 
of  avoiding  grossly  unphysical  results. 

The  conclusions  reached  regarding  the  two-point,  two-time  quasi-normality 
hypothesis  do  not  in  thfsnsolves  indicate  the  impossibility  that  the  simultaneous 
distribution  may  obey  a  quasi-normality  hypothesis.  Some  insight  into  this 
question  perhaps  is  provided  by  the  dissipationless  system  discussed  in  Section  5. 
For  this  system  it  was  seen  that  the  simultaneous  distribution  could  be  rigorously 
normal.  However,  the  third  moments  of  the  corresponding  two-time  distribution  do 
not  vanish,  nor  are  the  fourth  moments  of  the  two-time  distribution  normally 
related  to  second  moments.  This  may  be  seen  directly  by  expressing  the  derivatives 
33  ,  /3^,  3  T  -Xrt- /S^  at  'Z^  =  0  in  terms  of  the  simultsineous  distribution,  using 
the  equation  of  motion  (5»3),  and  verifying  that  they  do  not  all  vanish  for  a,[i,\,cr 
all  different.  If,  now,  viscosity  and  driving  forces  are  reintroduced,  the  third 
moments  of  the  simultaneous  distribution  take  on  non- vanishing  values.  It  seems 
at  least  plausible  that,  like  the  vanishing  of  the  third  moments,  the  normality 
relation  between  foiirth  and  second  moments  in  the  simultaneous  distribution  is  a 
property  peculiar  to  equi partition  and  may  not  be  realizable  when  there  is  a 
strong  mean  energy  transfer,  if  the  second,  third,  and  fourth  moments  are  to 
satisfy  the  conditions  for  the  existence  of  the  distribution  function. 

In  conclusion,  it  should  be  recalled  that  the  considerations  of  the  present 
paper  were  confined  to  examination  of  the  energy  balance  in  modes  for  which  the 
external  driving  forces  were  assumed  to  vanish.  The  results  obtained  do  not 
establish  the  impossibility  of  an  approximately  quasi-normal  distribution  for  modes 
in  the  energy  contsdning  range,  where  the  external  forces  cannot  be  ignored  in 
stationary  turbulence. 
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